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ABSTRACT 

Sagittarius A* in the Galactic center harbors a supermassive black hole and 
exhibits various active phenomena. Besides quiescent emission in radio and sub- 
millimeter radiation, flares in the near infrared (NIR) and X-ray bands are ob- 
served to occur frequently. We study a time-dependent model of the flares, 
assuming that the emission is from a blob ejected from the central object. Elec- 
trons obeying a power law with the exponential cutoff are assumed to be injected 
in the blob for a limited time interval. The flare data of 2007 April 4 were used 
to determine the values of model parameters. The spectral energy distribution 
of flare emission is explained by nonthermal synchrotron radiation in the NIR 
and X-ray bands. The model light curves suggest that electron acceleration is 
still underway during the rising phase of the flares. GeV 7-rays are also emitted 
by synchrotron self-Compton scattering, although its luminosity is not strictly 
constrained by the current model. If the GeV emission is faint, the plasma blob 
is dominated by the magnetic energy density over the electron kinetic energy 
density. Observations in the GeV band will clarify the origin of the blob. 

Subject headings: acceleration of particles - - black hole physics - - Galaxy: 
center — radiation mechanisms: nonthermal 
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Introduction 



Various observations have confirmed that the Galactic center, Sagittarius A* (Sgr A*), 



contai ns a super massive black hole of mass ~ 4xl0 6 M Q (e.g-. lGhez et al. Il2008r iGillessen et al. 
2009 aJJb). It has been found that Sgr A* emits radiation from radi o through X-rays (see 



Melia fe Falcke 



2001 



Melia 



20071. for review) and even TeV 7 -rays ( ITsuchiya et al. 1 12004 ; 



Kosack et al. Il2004t lAharonian et al. 1120041 : lAlbert et al. 1120061 ) . The bolometric luminosity 
of Sgr A* is ~ 10 36 erg s _1 and the emission is dominated b y radio. The spectral flux sharply 



falls off above v ~ 10 12 Hz ("the submillimeter bump") flZvlka et al. Ill995l : iFalcke et al. 



19981 ). In the quiescent state, X-rays at low luminosities were also observed by Chandra, 



which obtained the X -ray luminosity of Lx ~ 2.4 x 10 erg s in the 2-10 keV band 



( IBaganoff et al. 1 120031 ) . The quiescent state of Sgr A* ca n be described by a radiatively in- 



efficie nt accretion flow (RIAF) model (lYuan et al. 1120031 ) or a jet model ( IFalcke &: Markoff 



20001 ) . These models ascribe the submillimeter bump to synchrotron radiation by thermal 
electrons in magnetic fields of ~ 30 G. The X-ray emission in the quiescen t state is explained 
either by bremsstra hlung or by inverse Comp ton scattering in both RIAF ( lYuan et al. 1120031 ) 
and the jet models dFalcke k Markoff I [200oh . 



Strong flares occ ur frequently i n the X-ray band ( IBaganoff et al. I l200lh and near- 
infrared (NIR) band (IGenzel et al. 1 120031) . The det ailed properties of the flares of NIR 
and X-rays are reviewed in iDodds-Eden et al. X-ray and NIR flares occu r simul- 



taneously, with no significant delay (jEckart et al. Il2004t lYusef-Zadeh et al. 1 120061 ) . X-ray 



flares are always accompan ied by NIR flares, alth ough NIR flares may occur without an 
associated X-ray flare (e.g., iHornstein et al. 1 120071 ). It should be noted, however, that this 



might be instrumental, because the X-ray background in the Galactic center is proportion- 
ally larger than the NIR background and, thus, only bright X-ray flares are detected. The 
duration of NIR flare in the L' b and (3.80 /xm) observed in 2007 April 4 was about 100 
minutes ( iDodds-Eden et al. II2009I ). Simultaneously, X MM-Newton observed the flare in the 
X-ray band and the X-ray flare lasted about 60 minutes (IPorquet et al. II2008I ). The //-band 
light curve of the 2007 fla re has sub-structural variations on a timescale of ~ 20 minutes. 
Dodds-Eden et al. I ( 120091 ) attribute the NIR substructure to the fluctuation of the magnetic 
field. 

Observations of polariza tion of the IR flares a r e broadly consistent w i th that the IR flares 
are s y nchrotron origin (e.g.. lEckart et al. I l2006at iMarrone et al. I l2008t IDodds-Eden et al. 



2009; Yusef-Zadeh et al. 20091). On the other hand, the emission mechanisms of the X 



ray flares are still debat ed; upscattering of submillimeter photons (IMarkoff et al. 



Eckart et al. 



Yusef-Zadeh et al. 20061). sy nchrotron self-Comp ton (SSC) scattering ( Yuan et al. 



2004 



2001 



2003 



2006al Jbl; ISabha et al. 1 1201 00 , synchrotron emission from high energy 
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electrons (jYuan et al. 1120031 . 12004 ; iDodds-Eden et al. 1120091 ) . or inverse Compton scattering 
of NIR photons by ~ 10 M eV electrons responsible for the quiescent radio- millimeter emission 



2005 



(|Yusef-Zadeh et al. 1120091) . An orbiting hot spot model is also proposed (IBroderick k, Loeb 



hood ). This model requires the synchrotron cooling time longer than the orbital period 



to explain the light curve of NIR. However, the cooling time of X-ray emitting electrons is 
much shorter than the NIR emitting electrons. This makes the hot spot model inappropriate 
to the X-ray flare emission. 



Dodds-Eden et al. I ( 120091 ) explored the emission mechanisms of the 2007 IR/X-ray flare 



by various synchrotron and inverse Compton emission models. According to their model, 
SSC models need a large magnetic field such as 6000 G and a very small size of the flare 
emission region, e.g., ~ 0.0013-Rs, where i?g is the Schwarzschild radius. The magnetic field 
is too large compared with the inferred value during the quiescent state, i.e., 10 - 30 G, and 
the size is too small to account for the flare timescale of 100 minutes. On the other hand, the 
flare is well explained by "powerlawcool" model, in which electrons wit h a power law with a 



cooling break emit IR and X-rays by synchrotron radiation. Recently, ISabha et al. I (12 010) 



explor ed SSC models that include parameter regions not considered by IDodds-Eden et al. 
( 120091 ). Their model also requires a smaller source component size to give a sufficient inverse 
Compton scatter ing efficiency, so that m ulti-components are needed to obtain a broad time 



profile of a flare. iMarrone et al. I (120081 ) also constructed a SSC model of X-ray emission for 
a flare observed in 2006 July 17. The source size of their model is ~ Rs and the strength of 
magnetic fields is 1.5 G. The source size and the magnetic fields are smaller than those of 
our model by a factor ~ 10 as shown in Section [3J 

Because the origin of X-ray emission is still unclear as described above, in this paper we 
present a model that NIR and X-ray flares can be produced by synchrotron radiation from 
a single plasma blob. The simultaneity of NIR and X-ray flares favors a model of emission 
from a single region. The blob might be ejected from the region near the central black hole. 
As shown below, the size of the blob of our model is about 10 Rs and magnetic field is 
about 20 G. These values are consistent with the physical parameters of the accretion flow 
around the central black hole. In this model the flare duration is determined by the injection 
timescale of nonthermal electrons. The decay of the flares is owing to the radiative cooling 
and particle escape after turnoff of the injection of nonthermal electrons. In particular, we 
focus on the time-dependent behavior of the flare models. Because of various timescales such 
as cooling, injection, escape, etc., time-dependent simulations are necessary to investigate 
the flare properties in detail. It should be noted that most of the previous theoretical work 
on the flare emission has been restricted to the steady state models. 



Although we do not consider adiabatic expansion in this paper, there are phenomena 
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suggesting the effects of adiabatic expansion. For example, iHornstein et al. I ( 120071 ) found 
that their NIR measu reme nts are consistent wi t h a c onstant spectral index during a flare. 
Marrone et al. I (120081 ) and lYusef-Zadeh et al. I (120091 ). on the other hand, found that mil- 



limeter and submillimeter peaks delay the peak of NIR/X-ray flares. According to their 
models this delay is consistent with the adiabatic expansion of a self-absorbed source. The 
estimated expansion speed is much less than 0.1c. 

Our assumptions on the emission region are described in Section [21 Numerical results 
are given in Section El A summary of our results and discussion are given in Section [H 



Thro ughout this paper, we assume 8 kpc as a Galactic center distance (|Eisenhauer et al. 
2003) and 4 x 1O 6 M as the central black hole mass ( IGhez et al. 



2008; 



Gillessen et al. 



2009al Tbf). 



Blob Model 



We first assume that the emission region of the fl ares is a blob ejected from the central re- 
gion around the black ho le, e.g., from accretion flow (jWang. et al. Il2000t IWang fc Kusunose 



20021 lYuan et al. II2009I ). For simplicity of numerical calculations, we assume that the blob 
is a uniform sphere with radius R and magnetic field B. While RIAF or a jet supplies the 
quiescent component of the spectral energy distribution (SED) of radiation, mainly in radio, 
the temporal ejection of a blob produces flare emission. Since the blob is launched from the 
inner edge or near the inner edge of the accretion flow, the blob speed might be at least 
mildly relativistic. However, we assume that the blob moves at nonrelativistic speed, so that 
the relativistic effects such as beaming are neglected. Because relativistic jets have not been 
observed from the region near Sgr A*, it is possible to neglect the effects of the motion of 
the blob clS db first step of investigation. 

In our simulations the injection of nonthermal electrons in the blob triggers flares. Shock 
acceleration in the blob or magnetic reconnection can suppl y such nontherma l elect rons. 
We solve the kinetic e quations of electrons and photons as in Kusunose et al. I ( 120001 ) and 
Li fc Kusunose I (120001 ). We include the injection of nonthermal electrons, radiative cooling 
(synchrotron and inverse Compton scattering), and escape from the blob. The injection rate 
per unit volume and unit interval of 7, where 7 is the electron Lorentz factor, is given by 



9(7) 



0, 



exp( 



max ) > 



if 7 ^ 7min i 
if 7 > 7min , 



where K, p, 7 m i n , and 7 max are parameters. The value of 7 max can be a function of time, which 
allows different shapes of light curves. In our model described below, we assume that 7 max 
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increases linearly with time to fit the flare light curves of 2007 April 4. The normalization 
factor, K, is determined by assigning the value of the injection rate per unit volume, i.e., 



<?inj 



(2) 



7min 



A time-dependent model is also considered (see Section [3]). The injection duration of 
nonthermal electrons, ti n j, and the escape timescale of electrons from the blob, t esc , are also 
parameters. The emission mechanisms are synchrotron radiation by the nonthermal electrons 
and inverse Compton scattering of the synchrotron photons by the same electrons (SSC). 

In this paper we do not consider Comptonization of external soft photons such as radi- 
ation from accretion flow. We discuss the validity of this assumption in Section HI Further- 
more, the adiabatic expansion of an initially optically thick blob is excluded, because there 
is no significant de lay or asymmetry in the longer wavelength emission relative to the peak 
of the X-ray flare jDodds-Eden et al. lEoOflh . 



3. Numerical Results 



Numerical models are obtained for the parameters such as R, B, g in j, p, 7 m i n , 7 max , 
t; n j, and t esc - For a black hole mass Mbh = 4 x 10 6 M q , the Schwarzschild radius is Rs = 
2GMbh/c 2 = 1.2 x 10 12 cm. The light crossing time of this size is Rs/c ~ 40 s. The particle 
escape occurs by advection, so that we set t esc ~ several R/c. In our model the injection 
duration is strongly related to the flare timescale. The values of parameters are determined 
by the SED of IR and X-ray flares and light c urves . For this purpose, we use the d ata of the 
2007 IR/X-ray flare given in lPorquet et al. I J2008h and bodds-Eden et al. I J2009f ). 



3.1. Light Curves 

In our model, the emission in both IR and X-ray bands is produced by synchrotron 
radiation. The maximum value of vF v of synchrotron emission appears at ~ 10 16 Hz, which 
is in between the IR and X-rays bands, where F v is the energy flux per unit frequency. 
The luminosity of X-rays is sensitive to the values of 7 max and t in j. This is because X- 
rays are emitted by the electrons in the high energy end of the electron spectrum and the 
cooling time of those electrons is short. As soon as the injection rate of nonthermal electrons 
decreases, the X-ray luminosity declines. Because IR emitting electrons have longer cooling 
time than X-ray emitting electrons, the IR luminosity is more dependent on t e sc than the 
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X-ray luminosity. If t e sc is long enough, the X-ray emitting electrons begin to pile up as 
lower energy electrons that emit IR. 

We first examine light curves of the IR and X-ray bands for various models. For models 
in this subsection, we assume that R = 10 13 cm, B = 20 G, 7 min = 2, p = 1.3, t in j = 8R/c, 
and tesc = 5R/c. The blob is initially empty at t = and nonthermal electrons are injected 
for t > 0. In Figure [U we present light curves in the V (3.8 /im) and 2-10 keV bands 
for different models of 7 max and q irLy When 7 max is fixed, high-energy electrons are injected 
instantly at t = 0. Then the light curves rise steeply (dotted lines). It is found for this model 
that the light curve of X-rays increases faster than that of IR, which is not in agreement 
with the flare of 2007 April 4. To obtain slower rise of the light curves, the time dependent 
7max is assumed as shown by the dashed lines. Here we assumed that 7 max increases linearly 
with time from 500 to 5 x 10 4 during t — and t iny For t > t in j both models (the dotted 
and dashed lines) assume q^ = 0. By fitting model light curves with the data of the flares, 
we found that the above models have shorter decline timescales than the observed one. We 
then assume that the injection rate decreases gradually after t = t in j. As an example, we use 
a model with q(j) oc exp[— (t — t in j)/(^t in j)] for t > t iny In Figure [TJ a model with £ = 0.25 
is shown by the solid lines (model A). In this model 7 max is constant for t > t iny 

In Figures [2] and [3j we compare model light curves to the observed data for 2007 April 
4. Here the model light curves are the same as shown by the solid lines in Figure [TJ i.e., 
model A. The details of the model parameters are given in Table [TJ The size of the blob R 
is 10 13 cm and this corresponds to the light crossing time of 334 s. In this figure, the peak 
time of the luminosity of the model is shifted so that the model curves coincide with the 
observed light curves. We assumed the linear dependence of 7 max on time for t < t in ^ and 
the gradual decline of the injectio n rate for t > t\ ny The fac t that IR and X-rays attain the 



peak flux almost simultaneously (IDodds-Eden et al. 1 120091 ) implies that the value of 7, 



increases during the rising phase of the light curves, i.e., the acceleration of higher energy 
electrons is still underway in the early phase of the flare. 

The cooling time of X-ray emitting electrons is ~ 11 s and much shorter than the escape 
time t csc ~ 1.7 x 10 3 s, if 10 keV of X-rays are assumed. On the other hand, NIR is emitted 
by electrons with 7 ~ 10 3 . The synchrotron cooling time of those NIR emitting electrons is 
comparable with t esc , i.e., t coo i ~ 2 x 10 3 s. Thus the cooling is effective in X-rays during the 
decay of the flare and the decline time of the X-ray light curve is shorter than that of NIR. 
The decline of the NIR light curve is regulated by both t coo i and t esc . As shown in Figure [TJ 
a bump appears at 7 ~ 10 3 in the electron spectrum. This is because of balance between 
radiative cooling and escape. 
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3.2. Spectral Energy Distribution 

We next show the evolution of the flare SED of m odel A and c o mpar e the model to 



the observations in Figure HI The quiescent model by lYuan et al. I (120031 ) is also shown 
for comparison. The maximum flare luminosity is 5.0 x 10 37 erg s _1 at t — 9R/c. This 
luminosity is much higher than the submillimeter luminosity in the quiescent state, but this 
flare does not strongly affect the spectrum in the submillimeter wavelength. During the 
rising phase, synchrotron emission becomes luminous and attains the peak at t ~ 8R/c. The 
peak frequency increases because 7 max increases for t < t; n j. While the emission decays, the 
spectrum in the IR and X-ray bands become flatter and a dip appears at ~ 10 15 Hz. 

It is noted that the SSC component appears in the GeV region, whose intensity might 
be too small to be observed. However, other sets of parameter values allow a stronger GeV 
flux as shown below. 

Since the number of parameters is large and observed data are few, there is uncertainty 
in determining the parameter values to fit the data. We show another model (model B) with 
p = 1.8 in Figure [5j When p has a larger value, the value of 7 m i n should be larger for the 
radio flare not to exceed the quiescent radio emission. When the value of p is larger and the 
value of 7 m i n is fixed, there are more electrons with lower energy for the same IR and X-ray 
luminosities as the observed data. This constraint on 7 min may be not strong, if flares in 
radio occur simultaneously. 

As shown in Figure HI our model predicts the emission in the GeV band, which is 
produced by inverse Compton scattering of IR-optical photons, i.e., SSC process. The lu- 
minosity in the GeV band is not strongly constrained in our model. In Figure El we show a 
model with a larger luminosity in the GeV band (model C). A smaller value of the magnetic 
field and a higher injection rate produce a luminous SSC component. The difference of SEDs 
in Figures 0] and [6] also appears in the IR band of the decaying phase. In Figure 0] the IR 
spectrum becomes flat after t ~ 12R/c, while that in Figure E] keeps hard spectral shape. 



3.3. Electron Spectrum and Energetics 

The time evolution of the electron spectrum of model A is shown in Figure The 
number density of electrons increases for t < t in j because of injection. For t > t; n j, the injec- 
tion rate decreases exponentially with time. The higher energy electrons decrease because 
of radiative cooling. Also on timescale t esc , electrons escape from the blob. 

The time evolution of the energy densities of the radiation u ra d and the kinetic energy 
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of electrons u e is shown in Figure [8] for model A. In our model the magnetic field is assumed 
to be constant. The energy content is dominated by the magnetic energy density Ub, i-e., 
UB/u e ~ 16 at t ~ 8R/c. This is in contrast with blazar models, where the kinetic energ y 
of nonthermal electrons/positrons dominates over the magnetic energy ( IKino et al. 1 120021 ) . 
Blazars are relativistic jets aligned with the line of sight and the relativistic beaming effect 
enhances the luminosity, while our model of Sgr A* does not have relativistic motion. It 
should be also noted that blazars have a luminous SSC or external Compton component in 
the 7-ray band, while model A has only a weak SSC component. When the SSC component 
is more luminous as in model C, the electron energy dominates over the magnetic energy, 
i.e., ub/u £ ~ 0.1 at t ~ 8R/c. 



4. Summary and Discussion 

We calculated the SEDs and light curves from a blob as a model of the IR and X- 
ray flares observed from Sgr A*. Because we performed time-dependent calculations, it 
was possible to study the evolution of radiation and electron spectra in detail, which is 
different from most of previous work. Our model assumes that IR and X-rays are emitted 
by synchrotron radiation of nonthermal electrons injected in the blob. This model explains 
the characteristics of 2007 April 4 flare in IR/X-rays such as (1) the simultaneity of the 
X-ray and IR flares, (2) the steeper decline of the X-ray flare than the IR flare, (3) the flare 
timescale of 100 minutes. The flare timescale in our model is mainly determined by the 
injection time of nonthermal electrons, which is related to particle acceleration mechanisms 
and to be studied in future work. To compare the light curves we assumed that the value 
of 7 max increases linearly with time. If the value of 7 max is fixed during injection, our model 
does not explain the simultaneity of the X-ray and IR flares. The higher energy electrons 
responsible for X-ray emission are thus being accelerated in the rising phase of the flare. 

We compared the SEDs of different models. While model A has a faint SSC component 
in 7-rays and a larger value of ub/u £ , model C has a bright SSC one and a smaller value of 
us/u e . If model C is applicable, we will observe 7-rays in the GeV band. Future observations 
in the GeV band will constrain the model parameters as well as the energetics in the flare 
regions. Because the blob might be ejected from RIAF by magnetohydrodynamic processes 
and nonthermal electrons are accelerated in the collisonless shock, the observations of GeV 
7-rays will clarify the origin of the blobs and the particle acceleration mechanisms. 

The spectral index of NIR is almost constant during a flare in model C, while the 
spectrum becomes flatter during the decay in models A and B. In model C, the cooling time 
is longer than that of models A and B, because the magnetic field is weaker. Then only 
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escape time, which is independent of electron energy, regulates the decay of the NIR flare 
Previous work mainly modeled the SED of flares by steady state models a nd did not uti- 



l ize th e data of light curves. Time- dependent models were recently published by lSabha et al 



(120101 ) . They assumed SSC for X-ray flares and combined adiabatic expansion to account 
for light curves of 2007 April 4. To fit the data six source components are needed. While it 
is possible to attribute the IR substructure to those mult i- components of sources, it seems 
difficult to produce many emission regions coherently to cause flares. 

It is known that IR flares sometimes occur without X-ray flares. This might be owing 
to a small value of 7 max because of inefficient acceleration, rapid radiative cooling, or adia- 
batic cooling. If, for example, a blob is close to RIAF and external Compton scattering is 
important, synchrotron X-ray flares are suppressed. If this is the case, the effect of external 
Compton scattering appears as weak X-ray flares or hard spectral X-ray emission. This 
is because when synchrotron X-ray flares are absent, the kinetic energy of electrons with 
7 > 10 3 are given to the RIAF photons with v ~ 10 12 Hz, yielding X-rays with v > 10 18 Hz. 
The intensity of these X-rays depends on the distance of the blob from RIAF. 

We neglected the effect of the radiation from the RIAF in the flare emission spectra. By 
inverse Compton scattering of submillimeter photons with 10 12 Hz, electrons with 7 max ~ 10 4 
produces gamma-rays with ~ 1 MeV. If the submillimeter emission (L ~ 3.85 x 10 35 erg 
s _1 ) is received by the blob, the energy density of the submillimeter photons in the blob is 
~ 1.8 x 10 -2 erg cm -3 . Here we assumed that the photon density decreases as the square of 
the distance, that the distance of the blob from the submillimeter emission region is 1 AU, 
and that the submillimeter emission is isotropic. This energy density should be compared 
with that of synchrotron photons in the blob, ~ 0.3 erg cm -3 at t ~ 8R/c (model A). Thus 
external Compton scattering of submillimeter photons is negligible, if the distance is greater 
than about 1 AU. 



After completion of this work, we noticed a paper by iDodds-Eden et al. I (120101 ) . They 
solved a kinetic equation of electrons and obtained time-dependent flare spectra of syn- 
chrotron radiation: Note that they did not calculate SSC spectra. In their work flares are 
caused by the change in the injection rate and magnetic fields, while we assumed the change 
in the injection rate and 7 max - 
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Fig. I. — Light curves of the V band (3.8/zm) and the 2-10 keV band for various models. 
A model with a fixed value of 7 max = 5 x 10 4 is shown by dotted lines. The dashed lines are 
for a model with a time dependent 7 max . Here we assume that 7 max increases linearly from 
500 to 5 x 10 4 during t = and t- m j, where ti n j = 2.7 x 10 3 s. Both models assume that the 
injection is at a constant rate and stops at t = £ in j. A model with a time dependent injection 
is shown by solid lines. This model assumes that is constant for t < t in j and decreases 
as gi n j oc exp[— {t — ti n j) / (0.25t in j)] for t > t in j. The value of 7 max increases linearly as for the 
model shown by the dashed lines. For all these models, R = 10 13 cm, B = 20 G, 7 min = 2, 
p = 1.3, ti n j = SR/c, and t csc = 5R/c are used. 



2.0x10 




Fig. 2. — Light curves o f the V band. The flare data of the V band (filled squares) are 
from lDodds-Eden et al. I (120091 ). The solid line depicts model A. The value of 7 max increases 
linearly with time from 500 up to 5 x 10 4 . In Figures 2 and 3, the peak time of model light 
curves are shifted so that the peak time coincides with the observed peak time. 
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Fig. 3. — Light curves o f the 2-10 keV band. The flare data of the X-ray band (pluses) 
are from iPorquet et al. ] J2008h . Model A is shown by a solid line. 




Fig. 4. — Time evolution of the SED of model A. The SEDs are shown for t = R/c - 
19R/c at every 2R/c. SEDs in the rising phase are shown by thin lines and those in the 
decayi ng phase by thick lines. Radio to sub millimeter measurements are for the quiescent 



state jMarkoff et al. Il200ll:lzhao et al. ll20(iiUgp en circles). The flaring state in NIR (filled 



square) i s taken from iDodds-Eden et al. I (j2009h ■ The X-ray flare data (filled diamonds) 
are from iPorquet et al. I (120081 ) . TeV emissi on (lAharonian et al. 1 120041 ) is also shown by 



asterisks. The quiescent state RIAF model by hftian et al. I fhooal ) is shown by a dash- dash- 



dotted line for comparison. 
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Fig. 5. — Time evolution of the SED of model B with p = 1.8. Here, 7 m i n = 2 x 10 2 is 
assumed. The value of 7 max is increased from 1.3 x 10 3 up to 1.3 x 10 5 linearly with time 
during t = and i; n j. SEDs in the rising phase are shown by thin lines and those in the 
decaying phase by thick lines. 
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Fig. 6. — Time evolution of the SED with a luminous SSC component (model C). Here, 
B = 5 G and q m ^ = 7.5 are assumed. SEDs in the rising phase are shown by thin lines and 
those in the decaying phase by thick lines. 




Fig. 7. — Time evolution of the electron spectrum of model A. Here, n e (~f) is the electron 
number density per unit 7. The spectra are shown for every R/c. After t = 8R/c, the 
injection rate decreases exponentially. 




Fig. 8. — Time evolution of the energy densities of electrons (solid) and radiation (dashed). 
This model assumes that the magnetic energy density is constant, i.e., B 2 /(8n) = 15.9 erg 
cm" 3 with B = 20 G. 
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Table 1. Parameters 



Model B(G) p 7 min 7 max g in j(cm 3 s x ) 

A 20 1.3 2 5 x 10 4 1.6 

B 10 1.8 2 x 10 2 1.3 x 10 5 1.1 

C 5 1.3 2 5 x 10 4 7.5 



Note. - - All models assume R = 10 13 cm, t in j = 8R/c, and 
t eS c = 5R/c. The increase of 7 max during electron injection 
time is assumed and the values of 7 max listed above are those 



